Spiraling elliptic solitons in generic nonlocal nonlinear media 
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We have introduced a class of spiraling elliptic solitons in generic nonlocal nonlinear media. The spiraling 
elliptic solitons carry the orbital angular momentum. This class solitons are stable for any degree of nonlocality 
except for the local case when the response function of the material is Gaussian function. © 2013 Optical 
Society of America 
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Optical spatial solitons in nonlocal nonlinear media are 
attracting increasing attention during recent years in 
both theoretical [1-10] and experimental [11-14] aspects 
of research. The nonlocality plays an important role in 
the nonlinear evolution of waves. It may drastically mod- 
ify the properties of solitons. The solitons in bulk Kerr 
media may undergo catastrophic collapse [15,16]. The 
nonlocality of an arbitrary shape can eliminate collapse 
in all physical dimensions [6]. Nonlocality can support 
vortex solitons [17, 18]and multipolc solitons [19] which 
are unstable in local nonlinear media. 

In theoretical aspect, ellipse-shaped solitons have been 
reported in saturable nonlinear media, such as elliptic in- 
coherent solitons [20], elliptic dark solitons [21], and spi- 
raling elliptic solitons [22]. In experimental aspect, co- 
herent elliptic solitons [23] in lead glass which is nonlo- 
cal nonlinear media and elliptic incoherent spatial soli- 
tons [24] in photorefractive sceening nonlinear media are 
observed. 

In this Letter, we use the variational approach to de- 
rive the analytical spiraling elliptic solitons solution in 
generic nonlocal nonlinear media. We analyze the poten- 
tial function to study the stability properties of the class 
of solitons. 

The propagation of the optical beams in the nonlo- 
cal cubic nonlinear media can be modeled by the follow- 
ing generic dimcnsionless nonlocal nonlinear Schrodinger 
equation(NNLSE) [3,7], 
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where i/' = il>[x,y,z) is a paraxial beam, z is the 
longitudinal coordinate, V'^ = + '^v' ^ ^^'^ V ^'^'^ 
the transverse coordinates. An = J J ~ ^' iV ~ 
y')\i}}{x' ,y' , z)\dx' dy' is the normalized nonlinear per- 
turbation of refraction index, and R is the nonlinear 
response of the medium which is normalized, real and 
symmetric such that J J R(x,y)dxdy = 1. We suppose 
the material response to be Gaussian function [7,25], 
i.e.R{x,y) = l/{2n'wl,)exp[-{x^ +y^)/2wl^], where w„i 
is the normalized characteristic length of the material 
response function. 



By the variational approach [26], Eq.(l) can be inter- 
preted as an Euler-Lagrange equation corresponding to 
a vanishing variation 



l(-ip,il>*,'4>;,,ipl,ip^,-ip*,ipy,'ip*)dxdydz = 0, 



where the Lagrangian density / is given by [27, 28] 



i , ^^dip dip* I dip 2 , i^V'|2n 



(2) 
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R{x^^,y-7^)m,r^)\'-d^dTi. (3) 

We introduce a trial function, 

^{x, y, z) - A{z)G[X/b{z)]G[Y/c{z)] exp(z0), (4) 

where the Gaussian envelope is G{t) = exp(— 1^/2) the 
phase is = B{z)X'^ + Q{z)XY + Q{z)Y'^ + ip{z), and 
X = xcos/3(z) -t-2/sin/3(z),y = — xsin/3(z) -t- ?/ cos /3(z). 
Corresponding to the trial function we can obtain its 
power, P = J J \il;{x,y)\'^dxdy = nA^bc and orbital an- 
gular momentum(OAM), M = Im / ^ %p*{r x V-0)d^r = 
l/2F(fe^ — c^)8 with r = xe^ + yey Substituting the 
trial function above to the variational principle Eq.(2), 
we obtain the reduced variational equation 



Ldz = 0, 



(5) 



where L = J J Igdxdy, and Ig denotes the result of insert- 
ing the Gaussian ansatz (4) into the Lagrangian density 
(3). It also can be shown that the Hamiltonian corre- 
sponding to Eq.(l) is of the following form 



H 



dy' ' 2 
R{x-^,y-7j)\iPi^,r,)\^d^dTj 



dxdy. (6) 



After some algebraic calculations, L and H can be an- 
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alytically determined as 



L = 



46c 



{b^ + wl) {c-' + wl) 



(7) 



H 



b^ + c^ + ib'^B^c^ + ib^c^Q^ + b^c^e^ 

• (8) 



A^b^c^y/ib"^ + w2J (c2 + wl^) 



ib^+wl){c^ + wl) 



Following the standard procedures of the variational 
approach [26], we have b' = 2bB,c' = 2cQ,P' = 
(&2 + c2)e/(&2-c2) ,P' = 0,H' = and M' = 0. 
Primes indicate derivatives with respect to the evolu- 
tion variable z. So we can rewrite the Hamiltonian of 
the system as follows, 



H 



P 



(6'2 + c'2 + n), 



(9) 



Fig. 2 and Fig. 3 We introduce a nonlocal parameter 
a = max^Wm/b, Wm/c) to define the degree of nonlocal- 
ity for the beam in nonlocal nonlinear media. The larger 
is the nonlocal parameter, the stronger is the degree of 
nonlocality. In Fig.l, Fig. 2 and Fig. 3 u>,„ is 15.0, and the 
degree of nonlocality a is 7.5. 

An important aspect of any family of soliton solutions 
is their stability properties. We can study the stabil- 
ity characteristics of our analytical soliton solution by 
means of the analysis of the potential function 11(6, c). So 
we search the second derivative of the potential 11(6, c) 
with respect to b and c, then substituting Eq.ll and 
Eq.l2 into it we get 



962 



2 (6^ + c^) 



64c4 (62 + wl,) [64 + Gb^c"^ + c^ + A{b^ + c^) w^] 
b^c^ (6^ + 1462c2 + c") + [b^ + ISb'^c' + iib'c^ 



dbdc 



-^40*^)^^ + (6^ -f 56^0^ + 16c'') 

2 (6^ + c") 

" 63c3 [64 + 662c2 + c4 + 4 (62 + c2) 

6* + 146^ c^ + c^ + 12 (6^ + c^) 



(13) 



(14) 



n 



l_ 1_ 46"a2 4c2a2 

62 + c2 + (62 - c2)^ ^ (62 - c2)^ 
P 



(10) 



with a = M/P = 1/2(62 _ ^2)9^ 

Solitons can be found as the extrema of the poten- 
tial n(6, c). Letting 911/96 = and 9n/9c = 0, we can 
obtain the critical power and the 0AM. 



Pc = 



2(62 + c2)'7r[(6' 



- w: 



J] 



3/2 



64c4 [64 + 662c2 + c4 + 4 (62 + c2) wl] 

{b'-c'y[b'c'+{b'+c')wi] 

454^4 [54 + 662c2 + c4 -I- 4 (62 -f c2) wU ' 



(11) 

(12) 



We can also obtain the rotation velocity uj = /3' ~ 
2 (6^ 4- c^) (T^ (6^ — c^)^ . When the input power and 
0AM are chosen arbitrarily the spiraling elliptic soli- 
tons can be found the semi-axis of which are deter- 
mined by Eq.ll and Eq.l2. One example is shown in 
Fig.l with Pc = 127272.4 and = 0.560949 (other 
parameters arc 6 = 2.0, c = 1.0,0 = 0.373966, u;,„ = 
15.0 and w = 0.623277). The isosurface of intensity 
of the spiraling soliton is obtained from our varia- 
tional solution. Comparing two half widths obtained 
from variational solution, Wx — Vb^cos^ujz + c2sin2a;2: 
and Wy = \J c^cos^ujz + 62sin2cjz, with the numerical 
results we find an excellent agreement as is shown in 



d^U _ 2 (6^ + c^) 

9c2 ^ 64c4 (c2 + wl) [64 + 662c2 + c4 + 4 (62 + c2) 



b^c^ (64 + Ub'^c^ + c4) + (46^ + 3364c2 + 186' 



2c4 



+c«)^ 



(464 



5b^c^ + 4c4) 



(15) 



we can easily 
and A = 



From Eq.l3, Eq.l4 and Eq.l5 
get d^n/db^ > 0, d^U/dc^ > 
(92n/962)(92n/9c2) - (92n/969c)2 > when Wm ^ 
(the degree of nonlocality a is not zero). So 6 and c of 
the spiraling elliptic solitons what we have got analyt- 
ically by the variational approach are corresponding to 
the minima of the potential 11(6, c). So the soliton solu- 
tions are stable for any degree of nonlocality except for 
the local case. But we should mention that when the de- 
gree of nonlocality decreases the low-intensity oscillating 
tails which indicate the appearance of dispersive waves 
radiated by the soliton will occur [22,29]. The radiative 
tails take a portion of radiated 0AM from the soliton, 
and the reduction of 0AM in the main soliton leads to 
the slow reduction of ellipticity of the transverse rotating 
profile [22]. This can be verified in Fig. 4 and Fig. 5. In 
Fig. 4 and Fig. 5 w,n is 8.0, and the degree of nonlocality 
a is 4.0. 

In conclusion, we have obtained spiraling elliptic soli- 
tons in generic nonlocal nonlinear media by use of the 
variational approach. We show that this class of solitons 
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are stable for any degree of nonlocality except for the lo- 
cal case. Because of the appearance of dispersive waves 
radiated by the soliton the ellipticity of the spiraling el- 
liptic solitons will reduce when the degree of nonlocality 
becomes lower. Our theoretical results have been con- 
firmed by direct numerical simulations of the NNLSE. 
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Fig. 1. (color online) Propagation dynamics of the spi- 
raling elliptic soliton in nonlocal nonlinear media. The 
isointensity plot is at the level /m/2 of the elliptic soliton 
with Im = 20256.03 where /,„ is max|-0p.The normal- 
ized characteristic length of the material response func- 
tion Wrn is 15, and the degree of nonlocality a is 7.5. 




Fig. 2. (color online) Evolution of the beam width of the 
spiraling elliptic soliton in the direction of x axis. Num- 
berically obtained half width Wx (black line) is compared 
to the variational result (red line). The parameters Wm 
and a are 15 and 7.5 respectively. 




Fig. 3. (color online)Same as Fig. 2 but with plot cor- 
responding to the beam width of the spiraling elliptic 
soliton in the direction of y axis. 




Fig. 4. (color online) Evolution of the beam width of the 
spiraling elliptic soliton in the direction of x axis. Num- 
berically obtained half width Wx (black line) is compared 
to the variational result (red line). The parameters Wm 
and a are 8.0 and 4.0 respectively. 




Fig. 5. (color online) Same as Fig. 4 but with plot cor- 
responding to the beam width of the spiraling elliptic 
soliton in the direction of y axis. 
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